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ABSTRACT 

We measure Doppler shifts of interstellar absorption lines in HST/ STIS spec¬ 
tra of individual star clusters in nearby UV-selected galaxies. Values for systemic 
velocities, which are needed to quantify outflow speeds, are taken from the lit¬ 
erature, and verified with stellar lines. We detect outflowing gas in eight of 17 
galaxies via low-ionization lines (e.g., C II, Si II, A1 II), which trace cold and/or 
warm gas. The starbursts in our sample are intermediate in luminosity (and 
mass) to dwarf galaxies and luminous infrared galaxies (LIRGs), and we confirm 
that their outflow speeds (ranging from —100 km s _1 to nearly —520 km s _1 with 
an accuracy of ~80 km s _1 ) are intermediate to those previously measured in 
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Archive at the Space Telescope Science Institute, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 5-26555. These observations are associated with 
program GO-9036. 
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dwarf starbursts and LIRGs. We do not detect the outflow in high-ionization 
lines (such as C IV or Si IV); higher quality data will be needed to empirically es¬ 
tablish how velocities vary with the ionization state of the outflow. We do verify 
that the low-ionization UV lines and optical Na I doublet give roughly consis¬ 
tent outflow velocities solidifying an important link between studies of galactic 
winds at low and high redshift. To obtain higher signal-to-noise, we create a 
local average composite spectrum, and compare it to the high -2 Lyman Break 
composite spectrum. Surprisingly, the low-ionization lines show similar outflow 
velocities in the two samples. We attribute this to a combination of weighting 
towards higher luminosities in the local composite, as well as both samples being 
on average brighter than the “turnover” luminosity in the v — SFR relation. 

Subject headings: galaxies: ISM - galaxies: kinematics and dynamics - galaxies: 
starburst - ISM: outflows 


1. Introduction 

Galaxy - galaxy encounters likely trigger a significant fraction of all galactic star for¬ 
mation (e.g., Kennicutt 1989; Heckman 1998). The strongest starburst events often produce 
galactic-scale winds, which transport metals and dust from the interstellar medium (ISM) 
into the galactic halo, and possibly beyond (e.g., Larson 1974; Dekel & Silk 1988, Lehnert 
& Heckman 1996). These outflows are thought to heat and enrich the intergalactic medium 
(IGM) and even the intracluster medium. Although the amount of matter permanently es¬ 
caping the galaxy is still under debate, the interplay between star formation and galactic 
winds clearly plays a fundamental role in the evolution of galaxies. While feedback from 
active galactic nuclei may be important in larger galaxies (Scannapieco & Oh 2004), the 
deposition of energy and dispersal of metals by supernovae appear to shape the properties 
of smaller galaxies. In this paper, we explore the relation between outflow properties and 
galaxy morphology, environment, and star formation history via ultraviolet spectroscopy of 
local star-forming galaxies. 

Understanding the large-scale movements of material from the ISM into the IGM re¬ 
quires a panchromatic picture because the temperature of the outflowing gas spans a broad 
temperature range. Winds were first recognized in emission lines from warm, photoionized 
gas. Recent studies have observed cold gas via LIV/optical resonance absorption lines (Heck¬ 
man et al. 2000, Rupke et al. 2002, Schwartz & Martin 2004, Martin 2005). The warm and 
cold gas is thought to be swept-up by a much hotter wind (and/or entrained in it). X-ray 
observations detect the portion of the hot wind around 0.8 keV and find that this hot gas 
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contains the metals from the supernova ejecta (Dahlem et al. 1998; Martin et al. 2002). 

Much recent work has focused on absorption line studies of winds. Using the starburst 
itself as a background continuum source, the Doppler shifts of lines determine absolutely 
whether the gas is falling in (redshifted) or expanding outward (blueshifted). Moreover, un¬ 
like emission lines, whose strength scales as density squared, the strength of absorption lines 
are directly proportional to the column density of the ion along the sightline. Measurements 
for both local and distant starbursts often reveal blueshifted, interstellar absorption lines. For 
example, Schwartz & Martin (2004) measure outflows in nearby dwarf starbursts at speeds 
around —30 km s _1 . Heckman et ah (2000) find an average outflow speed of —100 km s -1 for 
luminous infrared galaxies (LIRGs). Martin (2005) and Rupke et ah (2002) detect outflow 
speeds of up to —700 km s -1 for two samples of ultraluminous infrared galaxies (ULIRGs). 
When these data for local starbursts are combined, the maximum outflow speeds measured 
in low-ionization lines are seen to increase with the global star formation rate (Martin 2005). 
The outflow speeds measured in distant starbursts (e.g. Shapley et ah 2003; Pettini et ah 
2002; Swinbank et ah 2004), redshifts z ~ 3, appear to be consistent with this empirical 
velocity - SFR relation. However, this comparison implicitly assumes that Na I absorption 
(rest-frame optical) is tracing the same outflow structures as the low-ionization UV lines 
measured for the high-redshift galaxies. 

Studying the kinematics of ultraviolet, interstellar lines in nearby galaxies serves several 
purposes. It provides a direct comparison between velocities measured from Na I and UV 
lines. It also increases the number of intermediate-luminosity starbursts with measured 
outflow velocities, which is important because the velocities from just three dwarf starbursts 
(Schwartz & Martin 2004) provide a lot of the leverage in the v — SFR relation (Martin 
2005). Additionally, the rest-frarne ultraviolet (UV) spectrum includes a wealth of resonance 
lines that sample a broader range of ionization states than do optical resonance lines (e.g., 
Kinney et ah 1993; Heckman & Leitherer 1997; Heckman et ah 1998; Kunth et ah 1998; 
Shapley et ah 2003; Chandar et ah 2004; Vazquez et ah 2004). In principle, measurements 
of ultraviolet resonance lines can determine how the outflow speed varies among different 
temperature components of the wind. For example, Vazquez et ah 2004 suggest an increase 
in outflow velocity with ionization energy in a high-resolution spectrum of a local dwarf 
starburst, NGC 1705. 

The UV cannot be observed in nearby galaxies from ground-based telescopes, so the 
spectroscopic capabilities of the Hubble Space Telescope provide unique data. As part of GO 
Program #9036, we examine UV resonance lines in Space Telescope Imaging Spectrograph 
(STIS) data of 17 nearby UV-selected galaxies. The least blended interstellar lines, which 
best constrain the gas kinematics, have first ionization energies less than that of H I. Hence, 



4 


we obtain the most information about the velocity of cold (i.e., ~ 10 2 K) and/or warm 
(< 10 4 K) interstellar gas, where much of the hydrogen may be neutral. In principle, hotter 
gas (> 10 4 K) is traced by C IV, Si IV, and N V, but stellar winds superpose broad absorption 
in these lines, making any interstellar component difficult to disentangle (Robert et al. 1993, 
Heckman et al. 1998). 

In this paper, we measure interstellar gas kinematics toward star clusters in a sample 
of nearby UV-selected galaxies. In §2, we discuss the HST/STIS observations, the data 
reduction, absorption line diagnostics, and systemic velocities. In §3 we discuss the ISM 
kinematics, including the kinematics of low-ionization absorption lines in individual galax¬ 
ies. In §4, we examine the outflow trends across the sample, including the dependence of 
outflow velocity on ionization energy, rotation speed, and star formation rate. In §5 we 
present a composite spectrum for the sample, examine the high-ionization absorption lines, 
and compare the composite spectrum to the high-redshift Lyman Break Galaxy compos¬ 
ite spectrum. The final section summarizes our results. We use heliocentric velocities and 
H 0 = 75 km s -1 Mpc -1 throughout this paper. 


2. The Data 

2.1. Observations, Data Reduction, & Spectral Extraction 

We have obtained Space Telescope Imaging Spectrograph (STIS) long-slit far- and near- 
UV spectra for 17 nearby starburst galaxies. We target individual UV-bright nuclear or 
near-nuclear (within 1 kpc of the galaxy center) star-forming clusters in the galaxies to pro¬ 
vide a strong background continuum source, against which to see absorption. The galaxies 
were observed using the 52"x 0.2" longslit and the G140L (FUV) and G230L (NUV) grat¬ 
ings, projected onto the 25 arcsec 2 MAMA detector. This configuration provides continuous 
wavelength coverage from 1175 A to 3100 A. The STIS MAMA detectors have an average 
pixel scale of 0.584 (1.548) A pixcU 1 for the G140L (G230L) detector. The spectral resolu¬ 
tion (FWHM), according to the STIS instrument handbook, is 8 pixels, or 0.2", for a fully 
extended source. 

The galaxies in this sample were selected on the basis of high UV luminosity; this is a 
tracer of star-formation. The sample was selected to include a range of galaxy morphologies 
and metal content (7.7 < log(0/H)+12 < 9.2). See Table 1 for salient properties of the 
sample galaxies. The young, hot stars in the clusters provide a strong UV background 
source against which we can see absorption by the ISM. Since these clusters are the most 
UV luminous, almost all of them have ages of 3-5 Myr (Chandar et al. 2004). The individual 
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star clusters in these galaxies are the same ones discussed by Chandar et al. (2004); they 
analyzed the absorption present in the stellar photospheres and winds to characterize the 
stellar population in each galaxy. We concentrate on the interstellar contribution to the 
spectra. Our sample differs slightly from the Chandar et al. (2004) sample. They include 
NGC 3049, NGC 5253 and Tol 89, which we do not because the data were taken with a 
different slit and are not part of this program. 

We define dwarfs to be those galaxies with M R > —18 (following the convention of 
Thuan & Martin 1981, Marlowe et al. 1997, and numerous other authors), and separate our 
sample into two subsets in which to study outflow kinematics: dwarfs and disks. However, 
since our sample is UV-selected, and thereby selected for the presence of active star-forming 
regions, we expect to see mergers and interactions in our sample, and the disk/dwarf clas¬ 
sification scheme is less than perfect. To this end, NGC 5102 is anomalous, since it is a 
lenticular galaxy, classified as type SA0 pec (i.e., not a dwarf) but it has an absolute blue 
magnitude of-17.11, which would place it in the dwarf subsample. Mkn 33 and Tol 1924-416 
both have Mg < -18, yet they are classified in the literature as blue compact dwarfs. Because 
of this dichotomy in classification, these three galaxies are separated from the disk and dwarf 
subsamples in this paper. 

The raw data files were retrieved from the Hubble Space Telescope (HST) archive. They 
were then co-added and processed through the CALSTIS pipeline, which performs global 
detector linearity corrections, dark subtraction, flat-fielding, wavelength calibration, and 
conversion to absolute flux units. We identify the same clusters as Chandar et al. (2004), 
and extract the spectra using the X1D package in IRAF 1 . This technique extracts a spectrum 
and subtracts a background spectrum, usually removing most of the geocoronal Lycc and O I 
A1302 emission. The spectra were normalized by dividing each G140L spectrum by the mode 
of the flux between 1250 A and 1500 A. The G230L spectra were normalized by dividing the 
spectra by a third-order polynomial fit between 2000 A and 3000 A. Since we are interested 
in the kinematics of the interstellar absorption lines, we retain the highest possible spectral 
resolution and do not rebin the spectra. 

Multiple clusters were observed in every galaxy except NGC 5102 and NGC 6764. If 
measured velocities for clusters in the same galaxy are the same (i.e., a difference in velocity 
cannot be measured within la), then the individual cluster spectra are co-added to produce 
a single spectrum with higher signal-to-noise. This is the case in all galaxies with multiple 


1 IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by AURA, Inc., under cooperative agreement with the National Science 
Foundation. 
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clusters extracted, except for NGC 3310. When multiple cluster spectra were averaged, 
the instrumental resolution for the individual clusters were weighted by flux and averaged. 
Table 2 provides the number of clusters averaged (Col. 2), as well as the maximum offset of 
the clusters from the slit center in either direction (Col. 3). (The slit positions themselves 
are presented in Table 2 of Chandar et al. 2004.) Spatial variation of the outflows along 
multiple sightlines is further discussed in § 3.2. 


2.2. Spectral Resolution 

The spectral resolution for a point source observed with the 0.2" ST1S slit is ~0.88 A at 
1500 A; for a source that fills the slit, the resolution is 4.67 A. Most of the clusters in our 
sample are neither point sources nor do they fill the slit. We measured the cluster surface 
brightness profiles from the two-dimensional spectra by summing up 500 columns (12 / .'5) 
along the slit spatially. We can equate the resulting spatial scale in pixels to an instrumental 
line width in Angstroms using the scale of 0.584 A pixel -1 in the G140L spectra. Most 
clusters have an instrumental FWHM between 1.5 A and 3.5 A, with a median size of 
2.7 A (or 607 km s -1 at 1335 A). These values (column 7 in Table 2) give a reasonable 
estimate for the resolution of individual cluster spectra (i.e., the instrumental width). 

The instrumental width is subtracted in quadrature from the measured width; line 
widths for the C II A1335 line are presented in column 8 of Table 2. For the four galaxies 
(NGC 1741, NGC 5996, NGC 7552, and Tol 1924-416) where we detect foreground/Milky 
Way halo lines, the instrumental line width is consistent with the widths in Table 2. 

The wavelength calibration will be affected if the extracted cluster does not fall in 
the middle of the slit (along the dispersion axis). From the STIS handbook, typical wave¬ 
length calibrations are accurate to ~0.2 pixels, or ~20 km s _1 at 1500 A. If the cluster is 
shifted by 1/4 of the slit width, or 2 pixels, then the wavelength calibration will be “off” by 
~200 km s _1 at 1500 A. This shift is large enough that it is easily discernible from wavelength 
measurements of the stellar lines. Comparison with stellar lines in our spectra to literature 
values confirm that the spectra used in this work do not suffer from any obvious zero point 
offset in the wavelength scale. 


2.3. Systemic Velocities 

Once the spectra are extracted, the kinematics are examined. In order to measure the 
kinematics of the warm/cold gas in the ISM, we need to compare the systemic velocity of 
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each object with the velocity of the low-ionization lines. The values of v sys for each object 
(i.e., the barycentric velocity combined with the rotation velocity at the cluster position) 
are given in Table 2. Ideally, the systemic velocity is derived from stellar lines in a high- 
resolution spectrum. Keck 2 /HIRES (Vogt et al. 1994) optical spectra of He 2-10, NGC 4214, 
and NGC 4449 provide a systemic velocity from the Mg I b-band absorption triplet, which is 
present in the K-type giants and supergiants in the starburst region. When high-resolution 
spectra of the star clusters are not available, a systemic velocity is taken from the litera¬ 
ture/NED. Preference is given to CO and H I velocities over emission-line velocities since 
the CO and H I likely trace the gas in the colder molecular clouds in which the starburst 
occurs, while the emission lines could be coupled to an outflow or other non-stellar process. 
When the value of v sys in the literature is not in agreement with the stellar lines seen in the 
UV spectrum, the stellar velocity is used; this is the case for NGC 3125, NGC 7552, and 
Tol 1924-416. (Stellar lines used include, but are not limited to, C III A1427, Fe V A1431, 
S V A1502, and He II A1640.) 

2.4. Spectral Line Identification 

Table 3 presents the wavelengths and ionization energies of lines measured in the spectra. 
Severely blended lines are noted. In the G140L/FUV grating, the best lines to use to 
characterize the purely interstellar gas (i.e., the least blended and strongest lines) are Si II 
(1260 A), C II (1335 A), Si II (1527 A), and Al II (1671 A). Although Ly a (1216 A) is an 
important diagnostic line, our local spectra are contaminated by strong airglow emission from 
the earth’s atmosphere, which makes it very difficult to disentangle/analyze the underlying 
Ly alpha emission and absorption from the starburst galaxy itself. The O I (1302 A) and 
Si II (1304 A) lines are blended, and the line combination can be affected by geocoronal 
O I emission if the background/sky subtraction is imperfect; therefore we exercise caution 
when fitting this line. We follow the convention of Shapley et al. 2003, and use an average 
wavelength for O I/Si II of 1303.27 A for this line. The G230L/NUV spectrum contains fewer 
lines, which tend to be weak. Mg I A2853, is a blend of stellar and interstellar absorption, 
though it is only detectable in six galaxies. The Mg II doublet near 2800 A is largely 
interstellar, and is detected in nearly all the galaxies though it is often difficult to deblend 
the doublet lines. A multitude of blended Fe II lines are present from 2344 A to 2600 A. 

2 The HIRES data referenced herein were obtained at the W.M. Keck Observatory, which is operated 
as a scientific partnership among the California Institute of Technology, the University of California and 
the National Aeronautics and Space Administration. The Observatory was made possible by the generous 
financial support of the W.M. Keck Foundation. 



The lines are all fit with Gaussian profiles using the splot task in IRAF. This task does 
not allow us to constrain the line fits (e.g., all Si II lines should have the same velocity 
and width), however it is a reliable method for fitting spectra at this resolution. To check 
the accuracy of our fitting, three galaxy spectra were also modeled with the SPECFIT 
program (Kriss 1994), which does constrain line parameters, and also produced an output 
model spectrum. The results differed by <20%. The lines were therefore fit with splot ; 
all low-ionization interstellar line velocities were averaged to give an overall low-ionization 
interstellar outflow velocity. 


3. Results 

The results are presented in Table 2. The resulting difference between the systemic/stellar 
velocity and the velocity of the low-ionization lines measures the outflow velocity of the cool 
gas entrained in the galactic wind (if any). Twelve of the 14 galaxies with detectable interstel¬ 
lar absorption lines show blueshifts, and two do not. The spectra of NGC 4449, NGC 4861, 
and VII Zw 403 are too noisy to detect any absorption. While NGC 4214 shows an outflow 
based on other data (e.g., Schwartz & Martin 2004), our data do not resolve this. The other 
four galaxies (NGC 1741, NGC 3125, NGC 6764, and Tol 1924-416) are consistent with zero 
velocity (within 1-cr). 


3.1. Doppler Shifts 

A low-velocity Doppler shift is difficult to detect, due to the blending of lines and the 
modest resolution of the spectra. The accuracy with which the absorption lines can be fitted 
depends on the noise in the spectrum; on average, we estimate an error of ~10% of the 
spectral resolution, or about 80 km s' 1 in fitting the kinematic center of a spectral line. 
Figure 1 shows the relation between signal-to-noise ratio per Angstrom, (S/N)^, and the 
fitting error, Sv: 

(7u(km s' 1 ) = —1.8 x ( S/N )^ + 104 km s' 1 (1) 

We find that for a typical S/N of 13, the error in fitting is 81 km s _1 , which agrees with our 
estimate. For a cleaner spectrum with S/N of 30 (per pixel), the error is 50 km s _1 , and for 
a noisier spectrum with S/N of 8, the error is 90 km s _1 . 

The velocities of all detectable interstellar lines are measured and the average outflow 
velocity Av out fi ow = v sys — vu ne is presented in Table 2. C II A1335, Si II A1260, and A1 II 
A1671 are the only three unblended lines which are present in every spectrum. Si II A1527, 
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while blended with the high-ionization C IV doublet, is detected and measured in all spectra. 
The C II and Si II spectra are presented in Figure 2. 

We detect outflowing gas in eight of the 17 galaxies. The detection threshold of 1-cr 
may seem rather liberal; however, a threshold of 2-cr would decrease the detection level from 
eight galaxies to seven. The galaxy which does not fall into the 2 -a level is NGC 5102. 

A velocity shift of about 80 km s” 1 will be detected, particularly if it is present in 
multiple interstellar absorption lines. If there are multiple lines of a species (e.g., Si II 
A1260, A1527), we know those lines should have the same velocity and width. Moreover, 
if ions have similar ionization energies (see Table 3), the corresponding absorption lines 
should arise from regions of gas at similar temperatures, and therefore should be spatially 
coincident. The outflow velocities (and errors based on the signal-to-noise in the spectrum) 
are presented in column 6 of Table 2. The highest-velocity outflows are clearly blueshifted 
from the systemic velocity. Projection effects due to a non-zero inclination angle (column 6 
of Table 1) are not considered because we do not know the opening angle of the outflow. 


3.2. Spatial Variation of Outflows 

The spatial variation in outflow speed along the slit is interesting since some starbursts 
show a more obviously global outflow. For example, many ultraluminous infrared galaxies 
show a coherent outflow over kpc scales (Martin 2005). Systems without an outflow observed 
near one cluster might show outflows along other sightlines. 

At high resolution, two distinct clusters in the dwarf galaxy NGC 4214 (NGC 4214- 
1 and -2, separated by ~ 500 pc) show different kinematics measured by optical Na D 
absorption (Schwartz & Martin 2004). NGC 4214-2 shows a blueshift of just -23 km s -1 , 
whereas NGC 4214-1 shows no shift at all, demonstrating the small (yet detectable and 
significant) differences in outflows between star forming clusters in the same galaxy. While 
this does not show that the outflows are tracking local chimneys from the disk, similar outflow 
velocities from spatially separated clusters can demonstrate the scale of spatial coherence 
in an outflow. Systems without an outflow observed near one cluster might show slightly 
different outflow kinematics along other sightlines, although at this resolution we cannot 
distinguish ~ 20 km s _1 differences. 

NGC 4214-2 was not observed with STIS. When NGC 4214-1 is fit with SPECFIT, 
the interstellar lines (C II, A1 II, Si II) converge to a fit where they are blueshifted by 
-32 km s -1 ± 49 km s -1 . Since there is no outflow detected in the optical spectrum of 
NGC 4214-1, it is not surprising that the UV velocity measurement is consistent with no 
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outflow. While by no means conclusive, this marginally demonstrates that outflows can be 
localized, although within the errors the UV data of NGC 4214-1 does not rule ont a slow 
outflow (as in NGC 4214-2). 

Multiple sightlines could be investigated in He 2-10. It has two main starburst regions, 
A and B, which are separated by 350 pc, and shows a large outflow velocity (Mendez et al. 
1999). Johnson et al. (2000) found similar outflow velocities from regions A and B. Our 
spectra from region A, the central starburst which contains five bright clusters, are shown 
in Figure 3. The regions all show similar kinematics, so we combine the spectra to increase 
the signal-to-noise. This is the case with every galaxy showing multiple clusters with the 
exception of NGC 3310. 

NGC 3310 is the only galaxy which exhibits significant (> la) kinematic differences in 
clusters - regions A, B, and C have different outflow velocities, and are therefore examined 
separately. The kinematics of these regions are presented in Table 2. One of the larger 
and more complex galaxies we observe is NGC 3310. The STIS slit is positioned to observe 
three large clusters near the center of the galaxy - see Figure 4. The de-redshifted spectra 
of regions NGC 3310-A, -B, and -C are shown in Figure 5. The spectra of the interstellar 
lines (dashed/red) make it clear that not only is there a significant outflow of interstellar gas 
from the star-forming regions, but this outflow changes across the slit. Region A is ~850 pc 
from region C, which is ~100 pc from region B. The H I velocities are the same in Regions 
B and C, whereas Region A is more redshiftcd by 30 km s _1 (Mulder 1995); this is too small 
a velocity difference to be detected in these spectra, as discussed above. However, there are 
large (>100 km s~ x ) differences in the outflow speeds from these three regions, with Region 
A (C) being the slowest (fastest) outflow; the slowest outflow is still significantly larger (Ar? 
= -339 km s -1 ) than any other Doppler shift in this sample. This suggests we are observing 
either local outflows from different regions of the disk (rather than some global outflow), or 
local variations within a global outflow. The scale over which kinematics differ is surprisingly 
small - regions NGC 3310-B and -C are separated by only ~100 pc. 


3.3. Comparison to Previous Results 

Previous measurements of outflow velocities from UV absorption lines are available for 
only three galaxies in our sample (e.g., He 2-10: Chandar et al. 2003, Johnson et al. 2000; 
Mkn 33 & Mkn 36: Kunth et al. 1998). He 2-10 is discussed below. The outflow velocity 
obtained for Mkn 33 agrees well with the measurement of Kunth et al. (1998). Mkn 36, 
however, does not agree: Kunth et al. (1998) find a redshift of +40 km s -1 whereas we 
find a relatively fast outflow. This may be due to the availability of many more lines in our 
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spectrum; Kunth et al. (1998) were only able to measure two absorption lines, 0 I A1302 
and Si II A1304. 

Larger galaxies with a more powerful starburst region tend to show a larger and therefore 
detectable Doppler shift (see §4.1 and §4.3). For example, NGC 7552 is a type SBab spiral 
galaxy. The STIS slit was placed across the nucleus. Figure 6 shows the absorption spectrum 
of this galaxy; a blueshift of —316 km s -1 ± 50 km s -1 is seen in the three unblended 
absorption lines. This is marginally consistent with the observations of Heckman et al. 
(2000), where the cold neutral medium is observed via Na D absorption, and is measured to 
have a blueshift of —216 km s -1 ± 20 km s -1 . 

The five separate clusters in region A of He 2-10 have the same outflow velocities to 
within < 30 km s -1 , (the fainter clusters have lower S/N and hence larger uncertainties). The 
low-ionization lines are fit with the SPECFIT program; this allows us to constrain the model 
fit, forcing lines of the same species to have the same velocity and width. Using the systemic 
velocity for this region from Mg-b absorption in high-resolution optical spectra (Schwartz 
& Martin 2005, in prep.), we fix the systemic velocity to be at the same redshift as the 
stellar (Mg-b) optical absorption. We measure an average UV low-ionization outflow speed 
of —170 km s -1 ± 8 km s^ 1 from the central starburst (commonly called Region A). This 
velocity is lower than the —360 km s -1 outflow (from C II and Si II) reported by Johnson 
et al. (2000). Part of the discrepancy, about 30 km s -1 , is due to our systemic velocities; 
but most of the difference is likely due to the uncertain position of the clusters in their 
larger aperture. In our recently obtained optical spectrum of He 2-10, the interstellar Na I 
absorption shows an outflow with four separate wind components (Schwartz & Martin 2006, 
in prep.). The highest velocity component is at A v^a = —128 km s -1 , which is consistent 
with the low-ionization UV absorption measurements presented here. 


3.4. Line Widths 

It is difficult to analyze the (Gaussian) full-width at half-maximum (FWHM) the same 
way we analyze the line velocities. When centering a line, there is usually an obvious ab¬ 
sorption minimum, as seen in the C II spectra presented in this paper. Therefore, measuring 
the kinematic center of a line is generally straightforward. However, it is more complicated 
to accurately measure the FWHM since the multiplets are often blended with other species 
arising from different ionization energies and therefore possibly different gas phases. For 
example, Si II A1527 (ionization energy 8.15 eV) is blended with the C IV AA1548, 1551 lines 
(ionization energy 47.89 eV), which often shows an extremely broad P Cygni profile. For 
this reason, we use the FWHM of the unblended C II A1335 line (ionization energy 11.26 
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eV) to characterize the FWHM of the gas in the cold/warm phase of the ISM (see column 8 
of Table 2). All references to FWHM herein imply a deconvolved, deredshifted width. 

The deconvolved line widths (measured at the Gaussian FWHM) measured range from 
fewxlOO km s -1 to > 1000 km s^ 1 ; these are broader than the Na D line widths seen in 
dwarfs (Schwartz & Martin 2004), LIRGs (Heckman et al. 2000) or in ULIRGs (Rupke et 
al. 2002, Martin 2005). The thermal line width is a few km s _1 for a gas at 10 4 K. The lines 
we observe are significantly broader than the instrumental resolution, although we do not 
know this resolution precisely. The absorption lines in the STIS spectra are not resolved. 

The estimated maximum/terminal velocity we see scales roughly as Av out fiow+FW H M /2] 
see Figure 7. (The exact fit is 0.62 ± 0.17x FWHM(C II). The profile shapes are described 
by a Gaussian profile. Generally, the lines get broader as the outflow velocity increases in 
magnitude. This is consistent with the physical picture of gas clouds being injected into a 
flow at the systemic velocity and being accelerated, as seen by Heckman et ah (2000). Since 
the instrumental FWHM comes from the spatial size of the cluster (see § 2.2), even with the 
largest instrumental widths we still see a non-zero FWHM in the C II line. This suggests 
that perhaps we are seeing a conglomeration of filaments and shells of interstellar gas at a 
variety of velocities, as seen in simulations (Fujita et ah 2005, in prep.), rather than a single 
cloud or simple shell. 


4. Kinematics of the Interstellar Medium — Analysis 

In Figure 8, we plot the outflow velocity Av against the cluster age, UV slope, and 
metallicity. There are no obvious correlations among these parameters seen in the data, 
whether or not the samples are divided into sub-samples by galaxy type (i.e. dwarf or disk). 
Therefore, outflow velocities are likely more dependent on the potential well of the galaxy and 
the star formation rate than the properties of the stars in the starburst. Surface brightness, 
star formation rate, and/or cluster mass may be the best parameters to examine, since the 
outflow seems more dependent on the galactic-scale environment than the stellar population. 


4.1. Rotation Speed vs. Outflow Velocity 

An important query is whether or not the outflowing interstellar matter is moving fast 
enough to escape the gravitational potential well of the host galaxy. We can parameterize 
the dynamical mass and therefore the gravitational potential of the galaxies with the galactic 
rotation speed. The escape velocity ranges from \[2v ro t (minimum) up to r-u 3v rot (for gas 
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extending ~ 3 kpc from the nucleus in an isothermal halo out to 100 kpc); a typical dwarf 
halo escape velocity is 100 km s _1 (Martin 1998). The resulting terminal velocities for 
the interstellar lines in our spectra are consistent with outflow being accelerated up to the 
escape velocity, but generally with no material exceeding v esc . We plot the outflow velocity 
as a function of rotation speed in Figure 9. The galaxies are separated into subsamples as 
described in § 2.1; the bright nearby dwarf starburst galaxy NGC 1705 (Vazquez et ah 2004) 
is plotted as well. 

Using the rotation speed as a proxy for halo mass, the overall trend suggests that more 
massive halos produce faster outflows. A least-squares fit gives a slope of 1.2 ± 0.2. The 
terminal velocity ( v term ~ Av ± 0.5xFWHM) may exceed the escape velocity significantly 
in two dwarf galaxies in particular - Mkn 36 and Mkn 209. This tentatively suggests that 
entrained outflow material from the ISM is likely to escape from these smaller galaxies. 
This result agrees with the debated theoretical picture that dwarf galaxies may be more 
susceptible to mass loss via superwinds, simply because they have a smaller potential well 
(e.g., Larson 1974; Dekel & Silk 1986; Martin 1998; Ott, Walter & Brinks 2005). 

As an example, the rotation speed of NGC 7552 from H I measurements is 230 km s -1 ; 
we estimate the terminal velocity of the gas in the outflow is Au + 0.5xFWHM ~ 700 km s -1 . 
This indicates that using a reasonable estimate for the escape velocity of v esc r3v rot (for 
an isothermal halo at 100 kpc and gas extending ~ 3 kpc; Binney & Tremaine 1987), the 
escape velocity is v esc ~700 km s -1 , and it is not likely that the gas is moving fast enough 
to reach the escape velocity for the galactic potential. 


4.2. Kinematics and Ionization Energy 

An important objective of this work is to check the consistency of outflow speeds mea¬ 
sured from different lines. We expect ions with higher ionization energies to trace higher- 
temperature gas, so hotter gas exhibits a faster outflow velocity (Vazquez et al. 2004). This 
relationship agrees with the scenario of hotter, less dense gas breaking free from a starburst 
region and expanding outward at a higher rate than the colder, more dense gas from the 
ISM that is loaded into the flow (e.g., Schwartz & Martin 2004 and references therein). 

High-resolution results (Vazquez et al. 2004) show a direct relationship between ion¬ 
ization energy and outflow velocity. In Figure 10, we show the outflow velocities for Al II 
and Si II versus C II velocity. Weighted least-squares fits to the Al II vs. C II data give a 
slope of 0.61 ±0.13. The correlation is noisy, but the Al II-absorbing gas seems to be moving 
slower than the C II-absorbing gas. The Si II-absorbing gas is less conclusive, with a slope 
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of 0.77 ± 0.15. However, the overall trend of A vcn > A vsui > A vaiii is consistent with 
IE(C II) > IE(Si II) > IE(A1 II) (where IE is the energy required to remove one electron and 
create the ion in question, as listed in Table 3. 

Optical interstellar Na I absorption lines have been detected for three galaxies examined 
- NGC 4214, NGC 7552, and He 2-10, as discussed in § 3.3. The low-ionization UV lines for 
NGC 4214 (Fig. 11; Schwartz & Martin 2004) and NGC 7552 (Fig. 6; Heckman et al. 2000) 
yield velocities slightly faster than the Na I measurements, although the optical spectra 
constrain the outflow speed better. The UV lines provide good constraints for He 2-10, 
where the average of several low-ionization lines is roughly consistent that of the highest 
velocity component seen in the Na I spectrum. The slightly faster speeds seen from the 
low-ionization UV lines (as compared to Na I) are consistent with the average I.E. of the 
UV lines being greater than I.E. of Na I. These results show that A v^ai < Avcii in galaxies 
where both lines are measured; this is in agreement with the absorption from different ions 
occurring in physically distinct regions of the outflow. 


4.3. Star Formation Rate vs. Outflow Velocity 

Another important parameter to investigate in relation to outflow velocity is the star 
formation rate (SFR). In Figure 12, we plot outflow velocity versus SFR, as measured by 
Ho luminosity (column 11 of Table 1), using the following equation from Kennicutt (1989): 

SFR (M n yr" 1 ) =- Llia - . (2) 

v 0 y ’ 1.26 x 10 41 V ; 

Although these Ho measurements were not corrected for internal extinction, the implied 
SFRs were found to be roughly consistent with those derived from the measured UV lumi¬ 
nosities. In contrast, the SFRs estimated from the FIR IRAS fluxes in the 60- and 100-micron 
bands were generally significantly lower than the UV-Ha -derived SFRs. Hence we adopted 
the Ha estimates of the global SFR. 

A general trend is apparent: overall, galaxies with a higher SFR exhibit faster outflow 
speeds. A simple least-squares fit gives a slope of 0.51 ± 0.59, This is consistent with the 
results for Na I velocity versus SFR presented in Martin (2005), and shows that UV absorp¬ 
tion lines and optical resonance lines yield consistent results. However, defining this relation 
empirically would benefit from more and higher resolution spectra. Also, since the SFRs are 
quoted for the entire galaxy here, more insight could be gained from a comparison to surface 
photometry of the individual star-forming regions. The total SFR is much greater than that 
of the cluster(s) observed in the largest galaxies such as NGC 3310. 
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The galaxies in this sample are intermediate in luminosity to the dwarfs in the Schwartz 
& Martin (2004) sample and the LIRGs in the Heckman et al. (2000) and Rupke et al. 
(2004) samples. The v — SFR relation of Martin 2005 predicts that velocities will level off 
around 10 1 —10 2 M 0 yr _1 , which is around the SFR of LIRGs. Therefore, the fitted velocities 
should be near or slightly less than the outflow velocities in the LIRGs; in general, they are. 


5. A Local Star-forming Galaxy Composite Spectrum 

To improve S/N, and study weaker lines as well as the systematic properties of the 
local star-forming regions, we created a composite spectrum of all the individual spectra. 
To produce our composite spectrum, we first shift all individual G140L spectra into the rest 
frame using the systemic velocities in Table 2. We then normalize each of the 16 individual 
spectra to the mode between 1250 A and 1500 A. The spectra are averaged and rebinned 
to a dispersion of 1 A. The composite local starburst spectrum is presented, along with the 
composite LBG spectrum from Shapley et al. (2003), in Figure 13. 


5.1. Kinematics of the Local Composite Spectrum 

The low-ionization absorption lines are blueshifted in the local composite spectrum. 

While the lines are still blended clue to stellar emission/absorption, there is little absorption 
redward of the systemic velocity, while blueward there is strong absorption. We used splot 
to fit Gaussian profiles to the low-ionization interstellar absorption lines (Si II, C II, Al II, 

O I/Si II) and found an average low-ionization outflow speed of A v out fiow = —142 km s -1 ± 80 km s _1 

In the composite spectrum, the higher S/N allows us to look for absorption from more 
highly ionized ions. There are three high-ionization lines: Si IV, C IV, and N V. N V is not 
present in the majority of the spectra, and when detected it is blended with Lyman a ab¬ 
sorption and corrupted by the residuals from sky subtraction of Lyman a emission. There is 
no distinct interstellar absorption component in the Si IV or C IV lines. The local composite 
C IV and Si IV are accurately described by a Starburst99 model spectrum (1-100 M 0 , O.25Z 0 
metallicity, instantaneous burst at 5 Myr; Leitherer 1999); we cannot confidently discern a 
high-ionization interstellar component in the composite spectrum in any high-ionization line. 
However, note that there is a detection of blueshifted interstellar N V, Si IV and C IV in 
the high-resolution spectrum of NGC 1705 (Vazquez et al. 2004); better resolution spectra 
of the clusters in this paper would likely yield a detection of some interstellar gas in these 
higher ions. 



5.2. Comparison with the High -,2 Composite Spectrum 


We compare our composite spectrum to the composite spectrum of the LBGs presented 
by Shapley et al. (2003). The samples differ in the number of galaxies included in the com¬ 
posite; Shapley et al. average 811 galaxies, whereas we only include 17 galaxies. Figure 14 
shows a histogram of the UV luminosities of the two samples; the local sample (average log 
L uv— 38.82 erg s^ 1 ) does not overlap at all with the LBG sample (average log L uv— 40.99 
erg s -1 ). Moreover, there is an “aperture bias” between the local and high -,2 samples. In 
the local spectrum we disproportionately weight the light from the clusters (which represent 
the youngest stellar population), whereas the high -z composite includes spectral features 
from both star-forming clusters and a diffuse “field” component. Chandar et al. (2004, 
2005) discuss in detail the contribution of the more diffuse “field” population to composite 
spectra. 

Given the differences in the two samples, the similarities between the local composite 
spectrum and the LBG composite spectrum are striking. First, the slopes in the low- and 
high-z composite spectra are remarkably similar: roughly -1.35 and -1.5, respectively (slope 
f3 between 1240 A and 1600 A for F\ oc \P). Using Starburst99 data (Leitherer et al. 1999), 
we can look at how the UV slope in a starbnrst changes with time. Both instantaneous and 
continuous bursts show a constant slope as a function of time out to ~10' 2 years, so the 
similar slopes are not hard to produce. The low-ionization interstellar absorption lines also 
show similar average outflow velocities: A v ou tfi ow — —142 km s _1 ± 80 km s^ 1 in the local 
sample, and A v ou tfiow = —150 km s' 1 ± 60 km s -1 in the LBGs. 

The similarity in low-ionization outflow velocity is quite surprising since all of the galax¬ 
ies that comprise the high-redshift composite spectrum have higher UV luminosity than even 
the most luminous member of the sample presented here. However, the velocity in the local 
composite is weighted towards the higher signal-to-noise spectra, i.e. the more luminous 
clusters, which have intrinsically higher outflow velocities. Therefore the dwarfs, with their 
overall lower outflow velocities, have less of an overall effect on the composite than the 
brighter galaxies, whose higher luminosity and higher outflow speeds are closer to that of 
the LBGs. While the composite outflow velocities may be similar in low-ionization lines 
between the two samples, the velocity distributions are different. Moreover, the SFRs for 
the two samples are both at or above than the “turnover” in the v — SFR relation of Martin 
(2005). Thus, an increase in SFR above this level ( 10 1 - 10 2 M 0 yr -1 ) will not result in 

an increase in the outflow velocity. 

Further examination of the composite spectra reveals that the local lines are somewhat 
broader than the high-z lines. This can largely be explained by a difference in spectral res¬ 
olution (local composite resolution ~800 km s -1 , z ~ 3 composite resolution ~600 km s^ 1 ). 
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The lines in the local composite are also sampling a broader range in luminosity, including 
both galaxies with low outflow velocity and those with large outflow. The combination of 
these two effects produces broader lines in the local composite spectrum. 

In particular, the C IV line in the local composite spectrum is broader than the ab¬ 
sorption in the LBG composite spectrum, even after accounting for the spectral resolution. 
We interpret this chiefly as an effect of age and aperture (since we are looking at younger 
stellar populations with STIS, as described above). Both the local and high -z composite 
spectra show a more pronounced P Cygni-like emission component in C IV than in Si IV. It 
is hard to reproduce the LBG C IV shape from Starburst99 models or local starburst galaxy 
components, but either the “blue edge” or the emission portion can be fit with these stellar 
models (Chandar et al. 2003). Hence, higher resolution observations are also required to 
definitively detect the high redshift outflows in high ionization lines. For example, the high 
resolution spectrum of the lensed z — 2.7 Lyman break galaxy MS 1512-cB58 (Pettini et al. 
2002) shows blueshifted interstellar N V, Si IV and C IV lines. 

As another interesting similarity between the local and high- 2 ' composites, some weak 
lines which appear to be noise in one of the spectra also are present in the other. The nature 
of these lines (e.g. the small lines at 1467 A, 1579 A, or 1689 A) is unknown, although they 
may be weak metal lines. Additionally, although the resolution of the composite spectra is 
too low to study column densities or covering factors via specific lines, the high resolution 
spectra of the local dwarf starburst NGC 1705 (Vazquez et al. 2004) and the Lyman break 
MS 1512-cB58 show similarly saturated absorption in the rest-frame UV. 


6. Conclusions 

We have detected many outflows of cool/warm interstellar gas in HST/STIS ultraviolet 
spectra of star-forming regions in nearby UV-selected galaxies. While the spectral resolu¬ 
tion of the G140L and G230L gratings is modest, eight of 17 galaxies show unambiguous 
blueshifted absorption in low-ionization interstellar lines at velocities from -100 km s^ 1 to 
-520 km s _1 . Somewhat lower velocities are measured for three of the dwarf galaxies, but 
higher resolution spectra are needed to confirm these estimates as typical measurement un¬ 
certainties are ±80 km s -1 . More massive galaxies with higher star formation rates generally 
host the higher velocity outflows. No redshifted interstellar absorption lines are found, which 
indicates a lack of inflows. That the lines are relatively broad (average C II line FWHM 
~ 800 km s^ 1 ) is likely indicative of sightlines probing a multitude of filaments and shells as 
seen in numerical simulations, rather than a single shell or bubble. 
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In practice, the velocities presented in this paper were heavily weighted by the strong, 
unblended lines from C II, Si II, and A1II. Since C and Si have ionization energies between 
those of Na and H, it is possible that (like Na I) C II and Si II are tracing gas where a 
significant fraction of the H gas is neutral. We looked for velocity differences among the 
UV lines with mixed results. The A1 II line tends to present lower velocities than the C II 
line. Since the ionization energy of A1 is less than that of C, this result supports previous 
suggestions that higher ionization gas moves faster than less ionized material (Vazquez et 
al. 2004). The ionization energy of silicon lies between that of aluminum and carbon, but 
the Si II line measurements show higher velocities than the C II line towards galaxies where 
Al II presents lower outflow velocities. Unfortunately, the high ionization lines from C IV, 
Si IV, and N V are too blended with stellar lines or airglow to reliably measure an interstellar 
component in individual spectra (but see the composite spectrum below). We conclude that 
the velocities will need to be measured more precisely before any definitive statements about 
velocity and ionization state can be made. 

To quantify spatial variations in the outflows, we measured line velocities toward all 
clusters intersected by the STIS longslit. For example, in NGC 4214 the Na I spectra clearly 
show outflow toward cluster ^2 but none toward cluster #1 only 500 pc away (Schwartz 
& Martin 2004). We found substantial differences in outflow velocities toward 3 clusters 
in an outer spiral arm of NGC 3310 on scales of 100 pc and 900 pc. In contrast, large 
velocities are measured all across the dwarf starburst He 2-10. Indeed, in all other galaxies 
where multiple clusters were observed, there is no observed spatial variation in outflow speed 
at this resolution. Although global winds, generally associated with nuclear starbursts, 
have received the most attention to date, our work draws attention to the need for more 
measurements of localized outflows in less extreme galaxies (Schwartz 2005, PhD thesis). 
Hence, the STIS slit may not have covered the star clusters with outflows in some galaxies. 
Considering also that dwarf galaxy outflows are generally slower than our detection threshold 
of 1-cr, the fraction of local star-forming galaxies with cool outflows is likely higher than our 
statistics (eight of 17) suggest. 

Studying galactic outflows in the ultraviolet facilitates direct comparison to the winds 
detected in high-redshift galaxies. The continuum slope of our composite spectrum for the 
local star-forming galaxies is indistinguishable from that of an analogous rest-frame UV 
spectrum for Lyman Break Galaxies (Shapley et al. 2003). The constituents of the latter 
sample are intrinsically brighter, but the two samples evidently have similarly young ages 
and average reddening. The low-ionization interstellar lines in these two composite spectra 
are remarkably similar with blueshifted absorption at about -150 km s -1 , which may reflect 
the flattening of outflow speeds seen as a function of SFR (Martin 2005). The lines for the 
local composite spectra are broader, reflecting a difference in spectral resolutions, as well as 
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the exclusion of the diffuse/held clusters in the local composite spectrum. 

The case for slower outflows in dwarfs rests on a very limited number of measurements. 
The STIS G140L spectra are not well-suited to kinematic measurements of dwarf starbursts. 
We are interested in the upper envelope of the velocity distribution over a broad range in star 
formation rates or, better yet, rotation velocities. Hence, the low velocities in NGC 4214 and 
NGC 4670 have little impact, as do the large uncertainties for NGC 5102. The main results 
are the revised outflow speed in He 2-10 (—170 ± 8 km s^ 1 ) and the new measurement for 
Mrk 33 (—200 ± 88 km s -1 ). These galaxies help fill in the gap between dwarf starburst 
and LIRG luminosities and are found to be consistent with the velocity limits described 
by Martin (2005). The large outflow speeds in Mrk 209 and Mrk 36 are of special interest 
since these galaxies are not very luminous. Better measurements of dynamical mass and star 
formation rate are needed for these galaxies. 
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Table 1. The Sample of Galaxies 


Galaxy 

(i) 

Galaxy 

Type 

(2) 

d 

(Mpc) 

(3) 

Vrot a 
(km s -1 ) 

(4) 

log O/H 
+12 
(5) 

i 

(deg) 

(6) 

M b 

(7) 

Age 

(Myr) 

(8) 

P 

(9) 

M(H I) 
(1O 7 M 0 ) 
(10) 

log 

L Ha 
(11) 

log 

L uv 
(12) 

References 

(13) 

He 2-10 

BCD 

9 

119±10 

8.9 


-17.32 

5+1 

-1.01 

30.5 

40.28 

42.74 

1, 2, -, 1, 3 

Mkn 33 

BCD 

19.5 

111+5 

8.4 


-18.31 

5+1 

-1.84 

62 

41.24 

42.61 

4, 5, 6, 7 

Mkn 36 

BCD 

6.9 

43+5 

7.8 


-14.02 

<1 

-1.83 

2 

39.85 

41.00 

4, 8, 4, 7 

Mkn 209 

BCD 

4.9 

47+21 

7.8 

30 

-13.30 

<1 

-2.45 

4.1 

39.90 

40.83 

9, 8, 10, 9, 7 

NGC 1741 

int. pair 

51 


8.2 


-21.15 

5+1 

-1.35 

1380 

41.77 

42.92 

-, 11, -, 13, 12 

NGC 3125 

BCD 

18.3 


8.3 


-17.81 

5+1 

0.62 


40.61 

42.02 

-, 14, -, -, 7 

NGC 3310-A b 

SABbc pec 

18 

311±10 

8.2 

52 

-20.13 

5+1 

-0.98 

588 

41.60 

42.73 

14, 15, 13, 14, 16 

NGC 3310-BC 1 

3 SABbc pec 

18 

311+10 

8.2 

52 

-20.13 

5+1 

-0.98 

588 

41.60 

42.73 

14, 15, 13, 14, 16 

NGC 4214 

IAB(s)m 

3.6 

70+10 

8.2 


-17.64 

4+1 

-1.20 

no 

>40.40 

41.23 

17, 1, -, 18, 5 

NGC 4449 

IBm 

3.6 


7.8 

51 

-17.79 

5+1 

-2.74 

135 

40.65 

41.19 

1. 17, 17, 18, 16 

NGC 4670 

BCD 

16 

110±21 

8.2 

28 

-17.93 

7+1 

-1.76 

110 

40.76 

42.31 

19, 13, 20, 20, 21 

NGC 4861 

BCD/Impec 

10.7 

54+3 

7.9 

82 

-17.25 



114 

41.09 

42.20 

5, 1, 6, 6, 6 

NGC 5102 

SA0 pec 

3.1 

95+12 

9.0 

70 

-17.11 

55-37 

1.50 

34 

>37.93 

40.72 

22, 23, 22, 22, 24 

NGC 5996 

SBc 

47 

142+10 

8.9 


-20.76 

4+1 

-1.29 

138 

41.16 

42.73 

25, 13, 26, 27 

NGC 6764 

SB(s)bc/Sy2 

32 

187±23 

8.7 

44 

-19.97 

3+1 

0.70 

380 


42.03 

28, 25, 28, 25, - 

NGC 7552 

SB(s)ab 

21 

230+30 

9.2 

31 

-20.36 

5+1 

-0.71 

870 

41.05 

42.29 

29, 13, 30, 31, 27 

Toll924-416 

BCD 

37 

45+10 

8.1 


-19.54 

<1 

-2.67 

300 

41.84 

42.95 

32, 13, -, 32, 6 

VII Zw 403 

BCD 

4.5 

20+3 

7.7 


-13.77 



6.9 

39.27 

40.46 

5, 8, -, 7, 7 


References. — (1) Sauvage et al. 1997, (2) Kobulnicky et al. 1999, (3) Johnson et al. 2000, (4) Bravo-Alfaro et al. 2004, (5) Davidge 1989, (6) Thuan 
et al. 2004, (7) Gil de Paz et al. 2003, (8) Perez-Montero & Diaz 2003, (9) Viallefond & Thuan 1983, (10) Stil & Israel 2002, (11) Guseva et al. 2000, 
(12) Dopita et al. 2002, (13) Heckman et al. 1998, (14) Mulder et al. 1995, (15) Pastoriza et al. 1993, (16) Kennicutt et al. 1989, (17) Martin 1998, (18) 

Karachentsev et al. 1999, (19) Hunter et al. 1996, (20) Hunter et al. 1994, (21) Marlowe et al. 1995, (22) van Woerden et al. 1993, (23) Chandar et al. 

2004, (24) McMillan et al. 1994, (25) Contini et al. 1997, (26) Kandalyan 2003, (27) Buat et al. 2002, (28) Wilcots et al. 2001, (29) Lehnert & Heckman 

1996, (30) Beck et al. 2002, (31) Claussen & Sahai 1992, (32) Ostlin et al. 1998 


Note. — Col. (2): Galaxy type is from RC3. Col. (3): Distance given assuming H 0 = 75 km s —1 Mpc -1 . Col. (4): Rotation speed is from the first 
reference in the last column. All rotation speeds are from H I measurements, and are corrected for inclination when available. Col. (5): 12 + log(0/H) is 
from the second reference in the last column. Col. (6): Inclination is from third reference in the last column. Col. (7): Absolute blue magnitude Mb is 
calculated from the RC3 magnitude Bj ,. Col. (8): H I mass is from fourth reference in the final column. Col. (9): Age is from Chandar et al. (2004). Age 
is based on the comparison of STIS FUV spectra with single burst STARBURST99 population synthesis models (Leitherer 1999). Col. (10): Power-law 
index of the UV continuum (F oc A^) in the STIS G140L spectra over the wavelength range 1240 - 1600 A after correcting for foreground Milky Way 



reddening. (3 is measured to an accuracy of ±0.01. (From Chandar et al. 2004.) Col. (11): Ha luminosity in ergs s -1 is from the fifth/last reference in 
the last column; these values are corrected for Galactic extinction and [Nil] (AA6548, 6583) contamination, but not for internal reddening. Col. (12): UV 
luminosity in ergs s —1 is from IUE flux at 1900 A, except He 2-10 and NGC 1741, where L jjy is at 1500 A, from Johnson et al. (2000) and Conti et 
al. (1996), respectively. To account for the UV spectral slope, (3 (Col. 9), we followed Bell & Kennicutt (2001). 

a Rotation speeds are from H I rotation curves for NGC 3310, NGC 4214, NGC 4670, NGC 5102, and NGC 7552. All other rotation speeds are from 
AV 50 , the H I velocity width at the 50% level. 

b NGC 3310-ABC are H II regions within NGC 3310. NGC 3310-A (see Figure 4) is equivalent to NGC 3310-1 in Chandar et al. (2004), while NGC 
3310-B and -C are part of NGC 3310-all in that paper. 



25 


Table 2. Absorption Line Properties 


Galaxy 

(i) 

No. of 
Clusters 
(2) 

Spatial 

Extent 

(3) 

Vs y s 

(km s -1 ) 

(4) 

Vs y s 

source a 

(5) 

A ^outflow 
(km s -1 ) 
(6) 

FWHMj„ s i 

(A) b 

(7) 

FWHM(C II) 
(km s -1 ) c 
(8) 

EW(C II) 
(A) 

0 ) 

Reference 

(10) 

He 2-10 

5 

+475, -172 

869+3 

Mg I 

-170+8 

2.5 

1136 

5.0 

1 

Mkn 33 

2 

-079, -471 

1465+10 

opt. lines 

-200+88 

3.8 

555 

2.4 

2 

Mkn 36 

6 

-076, -173 

646+5 

Hi 

-213+68 

2.6 

1203 

3.4 

3 

Mkn 209 

3 

+270, -272 

288+15 

CO 

-191+92 

1.6 

928 

1.7 

4 

NGC 1741 

4 

+277, -371 

3937+45 

opt. lines 

16+76 

3.2 

345 

2.1 

5 

NGC 3125 

1 

-079, -172 

964+73 

UV lines 

8+73 

2.5 

445 

1.1 

6 

NGC 3310-A d 

1 

-473 -476 

1040+12 

Hi 

-339+72 

2.9 

1088 

4.5 

7 

NGC 3310-B d 

1 

+875, +871 

1010+12 

Hi 

-414+82 

2.9 

1278 

5.1 

7 

NGC 3310-C d 

1 

+772, +770 

1010+12 

Hi 

-561+83 

1.9 

1382 

4.7 

7 

NGC 4214 

1 

-078, -1"3 

304+1 

Mg I 

-32+49 

3.3 

549 

2.1 

8 

NGC 4670 

2 

+275, -176 

1070+21 

Hi 

-65+73 

3.3 

527 

1.6 

9 

NGC 5102 

2 

-079, -172 

470+12 

Hi 

-98+78 

2.3 

859 

3.3 

10 

NGC 5996 

2 

-078, -279 

3297+32 

Hi 

-289+59 

2.7 

735 

2.7 

3 

NGC 6764 

1 

-078, -172 

2420+20 

CO 

-9+50 

2.6 

653 

3.1 

11 

NGC 7552 

2 

+177, +079 

1457+108 

UV lines 

-316+108 

3.5 

856 

5.6 

5 

Tol 1924-416 

2 

+476, -176 

2935+70 

UV lines 

-47+70 

2.3 

391 

1.1 

10 


References. —- (1) Schwartz Sz Martin 2005 (in prep), (2) Lequeux et al. 1995, (3) RC3, (4) Young et al. 1995, (5) Vacca Sz Conti 1992, 
(6) this paper, (7) Mulder et al. 1995, (8) Becker et al. 1995, (9) Hunter Sz Hoffman 1999, (10) van Woerden et al. 1993, (11) Eckart et al. 
1996 


Note. — The number of clusters given in Col. (2) is the number of individual cluster spectra extracted and averaged to produce the final 
spectrum of the region. The spatial extent given in Col. (3) is the maximum distance from the center of the slit to the clusters both above 
and below the center of the slit, in arcseconds. A v out fi ow = v — v sys given in Col. (6) is an average of all detectable (> 3cr) interstellar 
absorption lines. Since C II is the only unblended line detectable in every spectrum, its Gaussian FWHM is given in Col. (8), as it is believed 
to be representative of the typical absorption line. NGC 4449, NGC 4861, and VII Zw 403 were left out of this table due to very low (<1) 
S/N. 

a Source of v sys is given as one of the following: Mg I b-band (stellar) absorption lines, optical emission lines, H I 21-cm, CO lines, or UV 
absorption lines. 

b The instrumental line width is measured from the two-dimensional spectra, as described in §2.2. 

c The “actual” FWHM of C II is measured from the spectra and deconvolved in quadrature from the instrumental width, as described in 
§ 2 . 2 . 

d Systemic velocities for NGC 3310 are from H I velocity field maps from Mulder et al. (1995) to account for rotation, and show that the 
outflow velocities in different star-forming regions of NGC 3310 are certainly not identical. 
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Table 3. Absorption Lines of Interest 


Line A (lab) a 

I.E. b 

Astrophysical 

Notes 

(A) 

(eV) 

Source c 



Blends with Ly a absorption 
Blends with Ly a absorption 


cm 

1175.65 

24.38 

Stellar 

nV 

1238.80 

77.47 

Blend d 


1242.78 

77.47 

Blend d 

cm 

1247.38 

24.38 

Stellar 

Si II 

1260.42 

8.15 

ISM 

o I e 

1302.17 

0.00 

ISM 

Si II e 

1304.37 

8.15 

ISM 

ClI 

1334.53 

11.26 

ISM 

Si III 

1343.41 

16.35 

Stellar 

Si IV 

1393.75 

33.49 

Blend f 


1402.77 

33.49 

Blend f 

Si III 

1417.24 

16.35 

Stellar 

cm 

1426.80 

24.38 

Stellar 

FeV 

1430.57 

54.80 

Stellar 

sV 

1501.76 

47.22 

Stellar 

Si II 

1526.71 

8.15 

ISM 

cIV 

1548.19 

47.98 

Blend f 


1550.77 

47.98 

Blend f 

Fell 

1608.45 

7.87 

ISM 

He II 

1640.34 

24.59 

Stellar 

Al II 

1670.78 

5.99 

ISM 

Fell 

2344.21 

7.87 

ISM 

Fell 

2382.77 

7.87 

ISM 

Fell 

2586.65 

7.87 

ISM 

Fell 

2600.17 

7.87 

ISM 

Mg II 

2796.35 

7.65 

ISMS 


2803.53 

7.65 

ISMS 

Mg I 

2852.96 

0.00 

Blend 5 

Nal 

5889.95 

0.00 

ISM h 


5895.92 

0.00 

ISM h 


Blends with Si II A1250 
Blends with C III A1247 
Blends with Si II A1304 
Blends with O I A1302 


Often shows PCyg structure 
Often shows PCyg structure 
Blends with Si IV in PCyg profiles 


Blended with C IV 
Often shows PCyg structure 
Often shows PCyg structure 


a All lab (air) wavelengths are from Vazquez et al. 2004 except: C III A1427, Fe V 
A1431, and He II A1640 (from Morton 1991); Si II A1260 and Fe II A1608 (from 
Shapley et al. 2003). 

b The energy given as “I.E.” is the energy required to remove an electron from the 
parent ion to create the ion listed. For example, I.E.(A1 II) is the energy needed to 
remove an electron from neutral Al and create the Al ii ion. Values are from Cox 
2000. 

c Possible sources are predominantly stellar (“Stellar”), predominantly interstellar 
(“ISM”), or a blend (“Blend”). 

d N V AA 1238, 1242 is severely blended with Ly a absorption in nearly all galaxies, 
as is also the case with the Lyman Break Galaxies (Shapley et al. 2003). 

e O I and Si ii are blended together at our spectral resolution. We use the same 
average wavelength as Shapley et al. 2003: 1303.27 A. 







f Si IV and C IV are assumed to be blends of a P Cygni emission/absorption 
profile, absorption in stellar photospheres, and possibly some interstellar absorption 
(e.g., Heckman et al. 1998; Shapley et al. 2003; Chandar et al. 2004). 

«Mg II is assumed to be predominantly interstellar based on the findings of 
Vazquez et al. 2004 wherein no stellar components were observed for these ions 
in a high-resolution spectrum of a nearby dwarf starburst. Mg I is present as an 
interstellar line in NGC 1705 (Vazquez et al. 2004), but is also commonly seen in 
stellar UV spectra (Fanelli et al. 1992). Since the optical Mg I triplet (5172 A) is 
stellar as well, we assume this line as a blend. 

h Data for the optical resonance doublet, Na I, also called the Na D doublet, is 
provided for comparison. Schwartz & Martin (2004; 2005, in prep) analyze the Na 
D absorption in NGC 4214, NGC 4449 and He2-10, as discussed in the text. 



Fig. 1.— Uncertainty (Sv) in velocity (km s x ) as a function of median signal-to-noise for the 
spectra fit with SPECFIT: NGC 3310-ABC, NGC 4214, NGC 7552, and He 2-10. (Fitting a 
spectrum with SPLOT does not give a useful measurement of error.) The data are fit (dashed 
line) with a least-squares fit, giving a slope -1.8 km s” 1 ± 0.2 km s _1 per unit signal-to-noise. 
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Fig. 2.— (a) Normalized, rest-frame STIS G140L spectra of the C II absorption line. The 
rest wavelength of C II (1334.53 A) is marked by the dotted line. This is one of the only 
unblended interstellar lines which appears in all spectra. The top three rows shows C II in 
all the galaxies with outflows; the bottom row shows C II in the four galaxies where little 
or no outflow is detected, (b) The same plot, but for the Si II absorption line at 1260.42 A. 
The two lines plotted illustrate the variation in velocity between lines. For example, some 
galaxies (e.g., He 2-10) exhibit obviously blueshifted lines with similar profiles in both C II 
and Si II, whereas other galaxies (e.g., NGC 4670) show differences between the two lines. 
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Fig. 3.— (a) Normalized STIS G140L spectra of the low ionization absorption line profiles 
used to measure the warm/cold gas outflow velocity in He 2-10. Top: O I/Si II A1303 and 
C II A1335. Bottom: Si II A1263 and A1II A1671. The solid line plotted over the data is the 
fit from the SPECFIT program. The dotted lines represent the wavelengths corresponding 
to the systemic velocity for each of the absorption lines, which are labelled with the species 
and ionization state. The systemic velocity used is 869 km s -1 , as discussed in the text; 
this is consistent with the stellar lines seen in this spectrum. The interstellar absorption is 
measured to have A v ou tfi ow = -170 km s -1 . (b) Top: Same as (a), for G230L grating lines 
Fe II AA2344, 2374, 2383, Mg II AA2796, 2804, and Mg I A2853. Bottom: Keck/HIRES 
spectra (~7 km s^ 1 resolution) of Na D AA5890, 5896. The interstellar Mg I and Mg II lines 
are fit to have the same blueshift of -170 km s^ 1 as the optical absorption lines (Schwartz 
& Martin 2006, in prep.); the optical lines are statistically well-fit (y 2 = 0.75). The Fe II 
lines show smaller blueshifts, if any at all; these are blends of multiplet lines and nearby 
absorption in the Galactic halo and clouds. 



Fig. 4.— Archival HST/ WFPC2 F606W (MAST dataset #U2E67Q01T) image of the central 
region of NGC 3310, with position of the STIS slit overlaid. The star-forming regions A, B, 
and C are marked. Regions A and C are ~850 pc apart, and regions B and C are separated 
by roughly 100 pc. 
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Fig. 5.— STIS G140L spectra of the three bright star clusters (A, B, and C) observed 
in NGC 3310. The spectra have been de-redshifted. The solid lines (colored blue in the 
electronic edition) represent the wavelengths of stellar lines, and they are Si III A1417, C III 
A1427, and S V A1502. These are fit to check the value of the systemic velocity used. The 
dashed (red) lines are predominantly interstellar in origin; they are Si II A1260, O I/Si II 
A1303, C II A1335, Si II A1527, and A1II A1671. These are used to fit the velocity of the gas 
in the interstellar medium. The dotted (green) lines are the Si IV and C IV doublets, which 
are blends of the gas in stellar photospheres, winds, and possibly some interstellar gas. 
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Fig. 6.— Normalized STIS spectra of the low ionization absorption line profiles used to 
measure the interstellar gas outflow velocity in NGC 7552. The solid line plotted over the 
data is the fit from the SPECFIT program. The dotted lines represent the wavelengths 
corresponding to the systemic velocity for each of the absorption lines, which are labelled 
with the species and ionization state. The systemic velocity is 1457 km s^ 1 (from stellar 
absorption lines in the G140L spectrum); this is consistent with the stellar lines seen in this 
spectrum. The average outflow velocity of cold/warm gas as measured by low ionization 
lines such as C II, Si II, and A1II is 316 km s -1 . 



Fig. 7.— Average outflow velocity (An) for all interstellar absorption lines as a function of 
the FWHM of the C II line, which is often the only detectable, unblended, purely interstellar 
absorption line. The error bars are estimated from the relation between velocity fitting errors 
and signal-to-noise as discussed in 3.1. The solid line (black in the electronic edition) is a 
least-squares fit of the data which shows that FWHM ~ 2A v out fi ow (actual slope of fit = 
1.6 ± 0.50), as discussed in the text. 



Fig. 8.— Dependence of starburst region age (top), UV continuum slope fa (middle), and 
oxygen abundance log(0/H) + 12 (bottom) on outflow velocity. Age and fa are taken from 
models calculated in Chandar et al. (2004); NGC 5102 has an age of 55_37, and is not 
included in this graph for the sake of clarity. Oxygen abundances given as log(0/H) + 12 
are from references given in Table 1; abundance errors are not specified. The given outflow 
velocities are averages of all available, unblended interstellar absorption lines (as discussed 
in the text). Velocity uncertainties are estimated using the signal-to-noise of the individual 
spectra as discussed in Section 3.1, and are given in Table 2. The filled triangles (colored 
red in the electronic edition) represent dwarf galaxies; the (green) filled squares represent 
Mkn 33, NGC 5102, and Tol 1924-416, which aren’t classified as dwarfs or disks in this paper; 
(blue) open stars represent disk galaxies. 
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Fig. 9.— Outflow velocity (Av out fi ow = vu ne — v sys ; from an average of all available interstellar 
absorption lines) vs. galaxy rotation speed for the sources showing outflows. Disk galaxies 
(NGC 3310-ABC, NGC 5996, NGC 7552) are represented by stars (bine in the electronic 
edition); galaxies which are dwarfs by morphology but not by Mg (Mkn 33, NGC 5102, 
Tol 1924-416) are represented by a filled square (green). Dwarf galaxies are marked by filled 
triangles (red). NGC 1705 is marked by a filled circle (cyan; Vazquez et al. 2004). The error 
bars from the STIS data represent the velocity fitting error as discussed in §3.1. The dotted 
(red) diagonal lines show the relations v te rm = V^v ro t (minimum escape velocity) and 3v rot 
(escape velocity at 3 kpc in an isothermal halo extending to 100 kpc). The solid (black) line 
shows a linear least squares fit to the data; its slope is 1.2 ± 0.2. 
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Fig. 10.— C II outflow velocity versus A1II (top) and Si II (bottom) outflow velocities. The 
ionization energy of carbon is 11.26 eV, whereas the ionization energy of aluminum (silicon) 
is 5.99 eV (8.15 eV). The diagonal (dotted) lines represent vcn = vaiii and vcn = vsui- 
The least-squares fits of the lines give slopes of 0.77 ± 0.15 (0.61 ± 0.13) for A1II (Si II). 
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Fig. 11.— (a) Normalized STIS spectra of the low ionization absorption line profiles used 
to measure the warm/cold gas outflow velocity in NGC 4214, Region 1. The solid (red in 
the electronic edition) line plotted over the data is the fit from the SPECF1T program. The 
dotted lines represent the wavelengths corresponding to the systemic velocity for each of the 
absorption lines, which are labelled with the species and ionization state. The systemic veloc¬ 
ity is 304 km s -1 (Becker et al. 1995), which agrees with stellar lines seen in optical (Schwartz 
& Martin 2004) and UV (this paper). The measured line velocity is —32 ±49 km s -1 , which 
is consistent with the findings of Schwartz & Martin 2004 that no outflow is detected in 
optical resonance absorption lines in Region 1. (b) Normalized Na I Keck I/HIRES optical 
resonance absorption spectrum for NGC 4214-1 (Schwartz & Martin 2004). The doublet 
lines are marked D1 (A5889.95) and D2 (A5895.92). (c) Same as (b), for NGC 4214-2. Cold, 
interstellar gas is seen in absorption and is blueshifted from the systemic velocity by 23 
km s -1 . 
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Fig. 12.— Outflow velocity (A v out fiow = v ime — v sys ; from an average of all available in¬ 
terstellar absorption lines) vs. star formation rate as calculated from Ha luminosity using 
equation 2. The error bars represent the velocity fitting error as discussed in §3.1. Data 
point symbols and colors are the same as in Figure 9. The Ha luminosities have not been 
corrected for internal extinction, which would increase the luminosity (and therefore the 
SFR) by ~50%. The data point with the lowest SFR, NGC 5102, is a lower limit on SFR 
because the Ha luminosity is that of the nucleus only. The solid line is a least-squares fit to 
the data (excepting NGC 5102); the slope is 0.51 ± 0.59. 


Fig. 13.— (Top) Composite rest-frame UV spectrum of 811 Lyman Break Galaxies from 
Shapley et al. (2003); resolution ~2.8 A. (Bottom) Composite rest-frame UV spectrum of 
the 16 star-forming clusters in our sample; resolution ~2.9 A. The solid (colored blue in the 
electronic edition) lines represent the wavelengths of stellar lines, the dashed (red) lines are 
are predominantly interstellar in origin, and the dotted (green) lines are the high-ionization 
Si IV and C IV doublets, which are blends of the gas in stellar photospheres/winds, and 
possibly some interstellar gas. 
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Fig. 14.— Distribution of UV luminosities (L uv] ergs s -1 A -1 ) for the local sample (L uv < 
10 40 ; solid bold line - red in the electronic edition) and the LBG sample (log L uv> 10 40 ; 
Shapley et ah 2003; dashed line - blue in the electronic edition). 
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